Prawdopodobnie rozkład długości hospitalizacji nie zależy od charakteru udaru, co oznacza, że powtórny udar leczy się w ten sam sposób (stosując te same procedury)

zgony <- paste0("pZgon", c(7, 30, 90, 365))

wiek_plec_m <- dane[!(transport) & PLEC == "M" & lUdar < 3 & lUdar == kUdar, 
                      .(pZgon7 = 100*mean(zgon7),
                        pZgon30 = 100*mean(zgon30),
                        pZgon90 = 100*mean(zgon90),
                        pZgon365 = 100*mean(zgon365)
                        ),  
                      by = c("WIEK", "kUdar")
                      ] 

wiek_plec_k <- dane[!(transport) & PLEC == "K" & lUdar < 3 & lUdar == kUdar, 
                         .(pZgon7 = 100*mean(zgon7),
                        pZgon30 = 100*mean(zgon30),
                        pZgon90 = 100*mean(zgon90),
                        pZgon365 = 100*mean(zgon365)
                        ),   
                      by = c("WIEK", "kUdar")]

# plot_ly(wiek_plec_7_m, 
#         x=~WIEK, y=~pZgon7, 
#         color = ~kUdar, 
#         colors = kolory, 
#         type = "bar") %>% layout(barmode = "overlay", 
#                                  bargap = .1, 
#                                  hovermode = "x"
#                                  )
# 
# plot_ly(wiek_plec_7_k, 
#         x=~WIEK, 
#         y=~pZgon7, 
#         color = ~kUdar, 
#         colors = kolory, 
#         type = "bar") %>% layout(barmode = "overlay", 
#                                  bargap = .1, 
#                                  hovermode = "x"
#                                  )


subplot(lapply(zgony, function(x) plot_ly(wiek_plec_k, 
        x=~WIEK, 
        y=~get(x, wiek_plec_k), 
        color = ~kUdar, 
        colors = kolory, 
        showlegend = "if"(x == zgony[1], FALSE),
        type = "bar") %>% layout(annotations = list(x = 0.2 , 
                                                    y = 1.05, 
                                                    text = paste0(letters[which(zgony == x)], 
                                                                  ") śmiertelność ", 
                                                                  stri_extract(x, 
                                                                               regex = "[[:digit:]]*$"), 
                                                                  "-dniowa") , 
                                                    showarrow = FALSE, 
                                                    xref="paper", 
                                                    yref="paper"
                                                    ),
                                 title = "Śmiertelność kobiet w zależności od wieku",
                                 barmode = "overlay",
                                 bargap = .1, 
                                 hovermode = "x")
                                 ),  nrows = 2)
subplot(lapply(zgony, function(x) plot_ly(wiek_plec_m, 
        x=~WIEK, 
        y=~get(x, wiek_plec_m), 
        color = ~kUdar, 
        colors = kolory, 
        showlegend = "if"(x == zgony[1], FALSE),
        type = "bar") %>% layout(annotations = list(x = 0.2 , 
                                                    y = 1.05, 
                                                    text = paste0(letters[which(zgony == x)], 
                                                                  ") śmiertelność ", 
                                                                  stri_extract(x, 
                                                                               regex = "[[:digit:]]*$"), 
                                                                  "-dniowa") , 
                                                    showarrow = FALSE, 
                                                    xref="paper", 
                                                    yref="paper"
                                                    ),
                                 title = "Śmiertelność kobiet w zależności od wieku",
                                 barmode = "overlay",
                                 bargap = .1, 
                                 hovermode = "x")
                                 ),  nrows = 2, margin = .1)

Powtórny udar dotyczy węższej grupy wiekowej. W grupie najstarszych pacjentów drugi udar nie zdarza się ze względu na wysoką śmiertelność. Można spróbować też dla grup wiekowych (2 różne podziały)

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in `[.data.table`(popWiekBK_M, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 1 column 'lPacj': 18.961507 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in `[.data.table`(popWiekBK_M, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 2 column 'lPacj': 12.561175 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in `[.data.table`(popWiekBK_K, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 1 column 'lPacj': 18.961507 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in `[.data.table`(popWiekBK_K, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 2 column 'lPacj': 12.561175 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in `[.data.table`(popWiekEks_M, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 1 column 'lPacj': 33.072466 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in `[.data.table`(popWiekEks_M, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 2 column 'lPacj': 4.642784 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in `[.data.table`(popWiekEks_K, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 1 column 'lPacj': 33.072466 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in `[.data.table`(popWiekEks_K, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 2 column 'lPacj': 4.642784 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels

To jest do omówienia z ekspertami i z prof. Kamińskim.